arXiv:astro-ph/0211598vl 27 Nov 2002 


Mon. Not. R. Astron. Soc. 000, 1-^ (2002) Printed 2 February 2008 (MN style file v2.2) 


The Tully—Fisher relation of cluster spirals at z = 0.83*t 

Bo Milvang-Jensen,^f, Alfonso Aragon-Salamanca,^ George K. T. Hau,^ 

Inger J0rgensen,^ Jens Hjorth^ 

^School of Physics and Astronomy, University of Nottingham, University Park, Nottingham NG7 2RD, UK 
^European Southern Observatory, Casilla 19001, Vitacura, Santiago, Chile 
^Gemini Observatory, 670 N. A’ohoku Place, Hilo, HI 96720, USA 

'^Astronomical Observatory, University of Copenhagen, Juliane Maries Vej 30, 2100 Copenhagen 0, Denmark 


27—Nov—2002, accepted for publication in MNRAS 


ABSTRACT 

We present the rest-frame i?-band Tully-Fisher relation for a sample of 8 cluster spiral 
galaxies at z = 0.83 and 19 field spirals at z = 0.15-0.90 based on VLT spectroscopy 
and HST photometry. No strong difference is detected between the cluster and the 
held galaxies, but we hnd some evidence that the cluster spirals are ~ 0.5-1 mag 
brighter than the held ones at a hxed rotation velocity. Although only a ~ 1.5-2cr 
result, if conhrmed with larger samples this effect could be due to the cluster spirals 
experiencing a period of enhanced star formation while falling into the cluster. 

Key words: galaxies: spiral - galaxies: evolution - galaxies: clusters: individual: 
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1 INTRODUCTION 


Ground-based and HST observations indicate that the disk 
galaxy population in rich galaxy clusters has experienced re¬ 
markable evolution since 2 = 1. It has been argued that the 
increase with time of the SO fraction and the simultaneous 
decrease in the spiral fraction suggest that star-forming spi¬ 
rals fall into distant clusters at a much higher rate than in 
the local Universe, and that these spirals ultimately become 



is being removed from the disk a brief period of enhanced 
star formation could be expected. The strong evolution of 
the cluster spiral population contrasts with th e mild evo¬ 
lutio n observed in the field spirals to 2 


1 (Vogt 


1999 


2001, and references therein). To quantify the evolution of 
the cluster spirals, we are conducting a programme measur¬ 
ing the stellar and dynamical masses and M/L ratios for a 
sizeable sample of morphologically-classified disk galaxies in 
rich galaxy clusters at 0.2 < z < 0.9. We here present the first 


* Based on observations collected at the European Southern Ob¬ 
servatory, Chile (66.A—0376). 

t Based on observations with the NASA/ESA Hubble Space Tele¬ 
scope, obtained at the Space Telescope Science Institute, which 
is operated by AURA, Inc., under NASA contract NAS5-26555. 

J E-mail: ppxbm@nottingham.ac.uk 


results on the Tully-Fisher relation (Tully & Fisher 1977) 
of cluster galaxies at 2 = 0.83. 

We assume Hq = 75 kms“^Mpc“^ and qo = 0.05 in 
order to easily compare with Vogt and collaborators. 


2 THE DATA 
2.1 Sample selection 


We targeted the galaxy cluster MS1054—03 at 2 = 0.83 since 
it is a rich X-ray selected cluster at high redshift and has 
extensive HST imaging and 8-m spectroscopy available. We 
based our galaxy selection on the spectroscopi c and photo¬ 
metric catalogu es of known cluster mem bers of van Dokkum 
([1999|) (see also [van Dokkum et al.|[2000[ ). The spectroscopic 


catalogue was based on an 7-band selected sample with 
20.0 < I < 22.7 (corresponding to —22.3 < Mb < —19.6 
at this redshift). The catalogue gives Keck-based H5 and 
[OH] equivalent widths (EWs) and spectral types (Emission, 
Absorption or E-l-A). The photometric catalogue gives HST 
based magnitudes, colours and morphologies. The combined 
catalogue contains 87 galaxies, for which 74 have both pho¬ 
tometric and spectroscopic information. From this sample, 
we selected spiral galaxies as follows: 


(i) Galaxies having spiral morphology and Emission spec¬ 
tral type (EW([OH]) ^ 5A). We included galaxy 2011, 
which has spiral morphology, although it was originally clas¬ 
sified as merger/peculiar (M/P). Two galaxies were removed 
due to geometric constraints (1733 and 1888X). This gave 6 
galaxies, of which 5 yielded rotation velocities. 
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Field spirals: 




FOl, 0.15, +30° F02, 0.18, +48° F03, 0.22, +48° F04, 0.23, -17° F05, 0.25, -62° F06, 0.26, -37° F07, 0.26, -7° F08, 0.29, -51° F09, 0.32, +43° 


% 






% 



FIO, 0.32, +23° 

Fll, 0.33, +35° 

r 

F12, 0.33, +41° 

F13, 0.37, -41° 

F14, 0.43, +29° 

F15, 0.47, -9° 

F16, 0.47, +30° 

F17, 0.49, -5° 

F18, 0.55, -10° 


F19, 0.68, -7° F20, 0.69, -30° F21, 0.76, -30° F22, 0.90, -10° 

Cluster spirals, and the cluster M/P galaxy 1801: 


'I V 

1403, 0.81, +20“ 1896, 0.82, +37° 2130, 0.82, +1° COl, 0.83, +3° 1763, 0.84, +37° 2011, 0.84, +21° 1459, 0.85, -10° 661, 0.85, -30° 1801, 0.83, 0° 

Figure 1. HST+WFPC2 F814W images of the galaxies with rletente d emission. R e low e a,nh image is given galaxy ID, redshift and slit 
angle (ffsliti measured counter-clockwise from vertical). The IDs follow [van Dokkum ( |i99e| ) for the cluster galaxies, with the exception of 
COl (a cluster galaxy found by us). The VLT slits were aligned with the major axes of the galaxies. The images shown have been rotated 
to the mask position angle and are 4" on the side. The intensity scaling is linear, and the intensity cuts are the same for all the galaxies. 


(ii) Galaxies having spiral morphology and a less secnre 
Emission spectral type (EW([OII]) SA, but with larger 
errors). This gave 4 galaxies, but only 2 rotation velocities. 

(iii) Galaxies having spiral morphology and no listed 
spectral type. After removing one galaxy (1354) due to ge¬ 
ometry, this gave 2 Sa galaxies. However, we were unable to 
measure rotation velocities for them. 


Twelve spiral galaxies in the catalogue (50%) did not have 
Emission spectral type, and these were not selected. All 
galaxies were required to have an inclination of i > 30° 
(with 0° being face-on). Finally, we included an M/P clus¬ 
ter galaxy with strong [Oil] emission (1801) out of curiosity, 
since it did not conflict with any of the high priority targets. 
Thus, we started with a sample of 13 known cluster galaxies. 
To £11 the remaining space in our two spectroscopic masks, 
additional galaxies were selected. We morphologically clas¬ 
sified the galaxies in the WFPG2 mosa i c whi ch were not in 
the cluster catalogue of van Dokkum (1999) and then se¬ 
lected galaxies with spiral/disk morphology and i > 30°. 
We imposed magnitude and colour limits corresponding to 
those spanned by the selected cluster galaxies, i.e., F814W 
= 20.1-22.7 and (F606W-F814W) = 0.8-2.2. When there 
was a geometrical conflict, the brightest galaxy was chosen. 
In order to fill all gaps in the mask we sometimes relaxed 
the magnitude and colour limits. There were relatively few 
galaxies in the red end of our search window. This supple¬ 
mentary sample of 34 spirals yielded one additional cluster 
galaxy (GOl, F814W = 22.25) for which the rotation velocity 
was measured, and 19 z = 0.15-0.90 field galaxies with ro¬ 


tation velocities. The magnitude range for these was F814W 
= 19.1-23.2. The field galaxies observed provided an ideal 
comparison sample since it was observed under the same 
conditions as our cluster sample. A montage of WFPG2 im¬ 
ages of the galaxies with detected emission is shown in Fig. 


2.2 Spectroscopy 

Observations were carried out with the FORS2 instrument 
at the VLT (UT2) on 2001 Feb 23. Two spectroscopic masks 
were designed with position angles (PAs) at right angles to 
each other, in order to be able to cover all galaxy PAs. The 
slits were aligned with the major axes of the galaxies, tilting 
the slits with respect to the PA of the mask by an angle 6siit 
(cf. Fig. . The slits were 1" wide in the dispersion direction. 
Typical slit lengths were 11”. The spectral resolution was 
FWHM = 4.2 A, and the pixel size 1.075 A x 0.201". The 
median wavelength range was 5400-7600 A. Each mask was 
observed for 7 x 30 min. The seeing in the combined frames 
was 1.04” for mask 1, and 0.94” for mask 2. 

Details of the spectroscopic data reduction will be given 
in Milvang-Jensen (2002, PhD Thesis, in prep.). Briefly, af¬ 
ter bias subtraction, cosmic rays were removed using the 7 
frames taken for each mask. The geometric distortion was 
mapped and removed using the edges of the spectra in the 
‘dome’ fiats. Dome and sky flats were used to correct for 
pixel-to-pixel variations and slit profiles respectively. The 
science frames were then cut up into individual slit spectra, 
wavelength calibrated, and sky subtracted. 
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H F06, 2^0.26, H/3, 27 A x 6.0" (21kpc) 
Vrot sin 104^ j km/s, i = 78.4° 



Fll, 2 = 0.32, [OIII]5007, 30A x 6.9" (28 kpc) 

Vrot sin i — 75^2 km/s, i — 44.6° 

rd,.p=c=0.67l“;“="(2.8l“;;kpc), rd,phot= a.ll^'Jkpc 



1403, 2 = 0.81, [Oil], 33A X 6.4" (42 kpc) 

14ot sin i = 270j] km/s, i = 72.3° 

rd,.pec=0.84j:“ “«" (5.61“ = kpc), rd,phot = 7.7l“ J kpc 





2011, 2 = 0.84, [Oil], 23A X 4.3" (29 kpc) 

Ffot sin i = 84^g km/s, i = 46.6° 

rd,.pec=0.17l“ “="(l.ll“ Jkpc), rd,phot= l.B^^ Jkpc 


661, z = 0.86, [Oil], 24A x 5.1" (34kpc) 

I4ot sin j —7lj[2i kni/s, i — 33.4° 

rd,.pec = l.lll“ f2"(7.4l=;5 kpc), rd.phot = 3.ol“- J kpc 



Figure 2. Illustration of the emission line fitting. The first column shows HST+WFPC2 F814W images of the 6 example galaxies, 
rotated to have the slit along the y-axis. The following columns show 2D spectral images: observed, best-fit model and residual, with 
wavelength along the x-axis. The intensity cuts have been adjusted from galaxy to galaxy, except for the residual images. The galaxy ID, 
redshift, line ID and spectral image size are given on the figure, as well as the fitted values of the projected rotation velocity (Frot sin i) 
and the emission line exponential scale length (rj spec). The inclination (i) and the F814W photometric scale le ngth (rj p^ot) are also 
given. For each galaxy the height in arcsec of all 4 images is the same. We rejected the fit of galaxy F06, cf. Sect. |2.3|. 


Depending on redshift, different emission lines were vi¬ 
sually identified in the 2D spectra ([Oil]AA3726.1,3728.8, 
Hy, H/I, [OIII]A4959, [OIII]A5007 and Ha). 


2.3 Rotation velocities 


Projected rotation velocities were determined from the 2D 
emission lin e spectra nsing the s ynth et ic rot ation curve 
method of ^imard fc Pritchet (1998, 1999). A simple 
parametrized intrinsic rotation curve and an exponential 
emission intensity profile are assumed. Synthetic 2D emis¬ 
sion line spectra are then produced taking into account 
known observational parameters such as seeing, spectral res¬ 
olution and slit-width. The model parameters are: projected 
rotation velocity Kotsini; exponential scale length rd,spec; 
total line intensity 7; residual background level 6; and [Oil] 
doublet intensity ratio R. The main inputs are a continuum- 
subtracted postage stamp spectral image of the emission 
line, an estimat e of the background noise, and the inclina¬ 
tion i (cf. Sect. 2.4). The Metropolis algorithm (Metropolis 


et al. 


1953) is used to search the parameter space to get 


‘best fit’ values and 68% confidence intervals. 
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The fitting was done for all well-detected emission lines. 
Examples of emission lines, model spectra and residual im¬ 
ages are shown in Fig. Out of the 31 galaxies with well- 
detected emission lines, 3 field galaxies did not yield accept¬ 
able fits. One (F04) had a strange fine morphology: the line 
showed a tilt only on one side of the nucleus, whereas on the 
other the line was flat. Two (F02 and F06) had an intensity 
profile that looked more extended than exponential. F06 is 
shown in Fig. The galaxies with rejected fits have been 
excluded from our analysis. The M/P cluster galaxy (1801) 
was surprisingly well matched to the model. However, the 
emission was awmmetric, being strongest in the two bright 
knots (cf. Fig. y). We show this galaxy in Fig. but we do 
not include it in the analysis. 


Two different intrinsic rotation curves were tried: a 
step-wise flat rotation cu r ve an d the ‘Universal’ rotation 
curve of Persic & Salucci (1991). The latter is a function 
of the luminosity of the galaxy. The differences in the model 
images are small due to the effect of the seeing and the 
relatively coarse sampling. For the low luminosity galaxies 
the Universal rotation curve gives marginally higher rotation 
velocities that the fiat one. For the high luminosity galax- 
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Figure 3. (a) High redshift cluster and field _B-band Tully—Fisher relation. Frrnr b ars nn Mt 

plot symbols. The solid line is the local Tully-Fisher relation from [Pierce fc Tully ( |l992| ) , cf. Vogt (1999): Mg = — 7.48 log Viot — 3.10. The 
dashed lines mark the Scr limits (cr = 0.41 mag). Galaxy 1459 is off the scale of the plot at (M]j , log VJ-ot) = (—22.70,0.99). (b) Residuals 
from the local Tully-Fisher relation versus redshift. The dotted line is a fit to the field galaxies (except galaxy F03): ATF = (—1.6±0.3)5:. 


ies (i.e. the cluster galaxies and some of the field galaxies) 
there is little difference between the results from two ro¬ 
tation curves. In the analysis we will use the results from 
the Universal rotation curve, since this rotation curve has 
some physical motivation. For galaxies with more than one 
emission line in the observed wavelength range, the rotation 
velocities derived from the different lines agreed within the 
errors, and weighted averages were used in the analysis. 

The ratio rd,spec/rd,phot had a median value of 1.1 and 
a typical range of 0.7-1.7, in broad agreement with what 


is foun ^ locally ( Ryder & Dopita 1994). This indicates that 
the fitting procedure is producing sensible results. The fitted 
rd.spec was typically in the range 0.2-1.2”, and almost always 
greater than zero at the 3cr level. Spectral signal is detectable 
and thus contributes to the fit up to typically > 2 scale 
lengths. 

2.4 Photometry 

Photometry was carrie d out on the F814W a nd F606W 
HST+WFPC2 images (van Dokkum et al. 1999). The typi¬ 
cal exposure time for the combined HST images was 3.3 ksec. 
Photometric zero points were taken from the May 1997 
WFPC2 SYNPHO T up date. The SExtractor programme 
(Bertin & Arnouts 1996| ) was used to measure total mag¬ 
nitudes. Rest-frame B magnitudes were derived from the 
observed F814W and F606W photometry. F814W matches 
almost exactly the rest-frame R-band at 2 = 0.83, and thus 
the B magnitudes can be derived very accurately for galax¬ 
ies close to this 2 . F606W matches the B-band at 2 « 0.37. 
Thus, B magnitudes for all the galaxies in our sample can 
be interpolated with reasonably small uncertainties. The in¬ 
terpolation (and in a few cases, small extrapolation) was 
carried out using spectral energy distributions (SEDs) of 
local galaxies with different spectral/morphological types 


( [Coleman, Wu fc Weedmm 198C| ) . Eor each galaxy, the SED 
(or linear combination of SEDs) that yielded the observed 
(F606W—F814W) colours at the galaxy’s redshift was cho¬ 
sen for the interpolation. The formal uncertainty in this 
transformation is very small (rms ~ 0.01 mag), but that as¬ 
sumes perfect knowledge of the filter response in each band. 
Tests with different filters indicate that systematic uncer¬ 
tainties of ~ 0.1 mag could be present, but these are negli¬ 
gible in our analysis. Absolute B magnitudes corre cted for 
internal extinction. Mg, were calculated following Fully & 


Fouque (1985). The correction, which is a function of i, was 


in the range 0.32-0.96 mag, with a median value of 0.61 mag. 

The disk inclinations (i) were deter mined on the F814W 
images using the GIM2D programme (Simard et al. 2002). 
The uncertainties on i were included in the uncertainties on 
Vrot = {Viot sin i)/(sini). The inclinations were in the range 
33-82° (90° being edge-on), with a median value of 68°. 


3 DISCUSSION OF THE TULLY-FISHER 
RELATION 

The Tully-Fisher plot is shown in Fig. |(a). Galaxies 1459 
and F03 deviate strongly from the rest. Galaxy 1459 may 
be two galaxies or a chain galaxy (cf. Fig. |^. Galaxy F03 
has a very low surface brightness (the lowest in the sam¬ 
ple). These two galaxies will be excluded from the analysis. 
In Fig. H(a) we plot the l ocal Tully-Fish er relation (TFR) 


from Pierce & Tully (1992), as adapted by Vogt (1999). Most 
of our field galaxies and all our cluster galaxies fall on the 
high luminosity/low velocity side of this relation. The ab¬ 
solute magnitude residuals at fixed rotation velocity (ATF) 
are plotted against redshift in Fig. |(b). For the field galax¬ 
ies, ATF becomes more negative with 2 (94% significance 
from a Kendall’s tau test). Assuming that the errors in ATF 
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are the individual measurement errors plus 0.6 mag added in 
quadrature (a guessed intrinsic scatter in the B-band TFR 
for this redshift range), a chi-square fit for the field galax¬ 
ies gives ATF = (—1.6 ± 0.3)z with Xr ~ l-O- If taken at 
face value and interpreted as luminosity evolution, the ef¬ 
fects mag at z = 0.5, much larger than the ~ 0.2 mag 
Vogt (1999) found at similar redshifts. This interpretation 
assumes a non-evolving TFR slope, which is important since 
our field sample has a built-in positive correlation between 
luminosity and redshift. Our s ample does no t allow us to 
constrain the TFR slope, but Ziegler et al. (2002) found 
some evidence for a slope change at z ~ 0.5 for a sam¬ 
ple of 60 held spirals. Thus, the derived luminosity evo¬ 
lution of the held spirals must be regarded with extreme 
caution. One might worry that the low ATF values instead 
rehected underestimated rotation velocities due to not de¬ 
tecting the emission at sufhciently large galactocentric dis¬ 
tances. In that case low ATF values should be accompanied 
by low rd,spec/rd,phot values, and Kendall’s tan tests for the 
cluster and the held galaxies show this not to be the case. 


In order to compare the cluster and held spirals, we will 
use a constant TFR slope equal to the local value. We will 
make this comparison under two alternative (and extreme) 
hypotheses. First, we will assume that the zero point of the 
held TFR does not evolve with a and compare the residuals 
from the local TFR (ATF) for the held and cluster spirals. 
Second, we will assume that the zero point of the held TFR 
evolves with 2 as shown by the dotted line in Fig. ||(b), and 
compare the residuals from the local TFR corrected for this 
evolution, ATFcorr = ATF — (—1. 62 ). In both cases we com¬ 
pare the cluster sample {Ndus = 7) with the full held sample 
(A^fieid = 18) and a held subsample spanning the luminos¬ 
ity range of the cluster galaxies {Mg < —19.8; A^fieid = 7). 
Two statistical tests are carried out, a simple difference of 
the mean values (assuming rms/\/iV uncertainties) and a 
Kolmogorov-Smirnov (K-S) test. For comparison purposes, 
the K-S probability that the held and cluster samples are 
drawn from two diherent populations has been translated 
into a number of sigmas for a normal distribution. The re¬ 
sults of these comparisons are given in Table The diher- 
ence between cluster and held samples are larger and more 
signihcant using the hrst set of assumptions, but it is clear 
that the data suggest the cluster spirals are ~ 0.5-1 mag 
brighter than the held ones at a hxed rotation velocity. The 
only ot her published cluster TF study at intermediate 2 is 


Metevier, Koo fc Simarc| (2002), who found a larger 


that of 

TFR scatter for their 7 cluster spirals at 2 = 0.39, but no 
evidence for a zero point or slope change. The larger look- 
back-time of our study could explain the fact that we do 
hnd some luminosity evolution. 


The increment in luminosity that we have found could 
be the result of enhanced star formation in spiral galaxies 
falling onto the cluster. However, these suggestive results are 
far from conclusive since we have studied one single cluster 
and our current sample is small. It is clear that a similar 
study of a reasonable cluster sample, spanning a broad red¬ 
shift baseline, can provide strong empirical results on the 
evolution of the cluster spiral population. 
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Table 1. Tully-Fisher differences: cluster versus field 


Variable 

Sample 

(cluster) — (field) [mag] 

pK-S 

■^difF. distr. 

ATF 

All 

-1.32 ± 0.43 (3.1cr) 

96% (2.0cr) 

ATF 

High L 

-1.00 ±0.46 (2.20-) 

87% (1.50-) 

ATFcorr 

All 

-0.71 ±0.42 (1.7cr) 

83% (1.4cr) 

ATFcorr 

High L 

-0.67 ±0.43 (1.60-) 

87% (1.50-) 
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